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Abstract 

We present an i?-band image of FS Tauri, a 0''2-separated classical T Tauri bi- 
nary system, taken with the Coronagraphic Imager with Adaptive Optics (CIAO) on 
the Subaru Telescope. This new image, combined with Hubble Space Telescope / 
Advanced Camera for Surveys (H ST / ACS) F606W-band polarimetric images, shows 
that the binary has complicated circumbinary features, including a circumbinary disk, 
western and eastern arm-like structures, and two cavities at the northeast and south- 
west. The circumbinary disk is 630 AU in radius and the southeast side of the disk 
is bright in the i?-band. The brightness ratio (southeast/northwest) is 1.6 ±0.4. A 
single Rayleigh-like scattering model indicates that the disk is inclined by 30° to 
40° and that the southeast side corresponds to the near side along our line of sight. 
The if-band surface brightness of the southeast side decreases as r -L9±0,1 from 15.2 
mag arcsec -2 to 16.8 mag arcsec -2 . The outer portion of the disk is possibly more 
flared than its inner portion. The weak centro-symmetric polarization pattern and 
redder F606W— H color (4.2 ±0.2 mag) of the southeast side are probably caused 
by multiple scattering events from the dust grains associated with the binary. The 
F606W-band image shows the bright northwest side of the disk in contrast with the 
if -band image. The F606W— H color of the northwest side is between 1.7 mag and 
3.0 mag. We consider that Haro 6-5 B (FS Tauri B), 20" away, produces the neutral 
scattered light from the northwest side. This idea is supported by the polarization 
pattern of the northwest side, which is centro-symmetric with respect to Haro 6-5 B. 
The arms appear to encompass the western and eastern cavities, suggesting that the 
arms + cavity systems are created by a bipolar outflow from the binary. However, 
the direction of this outflow is inconsistent with that of outflows inferred from the 



circumbinary disk model. These differences may arise from misalignment between 
the circumbinary disk and the circumstellar disks. Another mechanism forming the 
arms + cavity systems is considered to be the inhomogeneous density distribution of 
materials in the circumbinary disk. 

Key words: star: individual (FS Taurus) — star: pre-main sequence — tech- 
niques: high angular resolution 

1. Introduction 

Many studies have showed that the majority of young stars (< 1 Myr) are associated 
with various circumstellar structures, such as envelopes, protoplanetary disks, and/or bipolar 
outflows. The disks are composed of gas and dust, the possible birthplaces of planets. The 
bipolar outflow creates a cavity in the envelope. The cavity allows optical and near-infrared 
(NIR) radiation to escape and illuminate an extended reflection nebula. These circumstellar 
features have been directly revealed using ground-based telescopes with adaptive optics (AO) 
and the Hubble Space Telescope (HST). 

A significant number of young stars appear in binary systems. NIR imaging surveys 
show multiplicity rates of 29 — 37% for Class I proto-stars in the p Ophiuchi molecular cloud 
(separation < 2000 AU; Haisch et al. 2002, 2004; Duchene et al. 2004). This rate is quite high, 
despite the relatively poor sensitivity at the small projected separations (< 110 AU). Some 
mechanisms are proposed for the formation of protobinaries, in which the core fragmentation 
process is important (e.g., Boss 2002). In the standard core fragmentation senario, a rotational 
and magnetic cloud core fragments into multiple dense objects and then a proto-binary is 
formed. The semi-major axis of the binary is determined by the initial core properties, especially 
by the magnetic and rotational energies and the angular momentum (Machida et al. 2008). 

In the more evolved stage (Class II or classical T Tauri stars), binary systems have two 
kinds of protoplanetary disks: circumstellar disks associated with each star and ring-shaped 
circumbinary disks around the binary systems (Artymowicz & Lubow 1994, 1996). The sizes 
of the circumbinary disks and circumstellar disks depend on the binary system's semimajor 
axis, eccentricity, and mass ratio (Artymowicz & Lubow 1996; Armitage et al. 1999). Several 
observations have detected scattered light from the circumbinary disks using optical and NIR 
imaging/polarimetric methods, for instance, around the GG Tau (Itoh et al. 2002; McCabe 
et al. 2002; Krist et al. 2005; Silber et al. 2000) and UY Aur (Close et al. 1998; Hioki et 
al. 2007; Potter et al. 2000) binary systems. Bipolar outflows also appear in several T Tauri 
binary systems (e.g., Hirth et al. 1997). However, observations have focused mainly on single T 

* Based on data collected at the Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 
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Tauri stars, although more than half of all T Tauri stars are binaries, according to observational 
results (separation < 1800 AU; Ghez et al. 1993; Leinert et al. 1993) and theoretical predictions 
(e.g., Nakamura & Li 2003; Machida et al. 2008). A limited number of studies have examined 
the structures around T Tauri binary systems. 

FS Tauri A (hereafter FS Tau or FS Tau binary) is a classical T Tauri binary system in 
the Taurus molecular cloud (d ~ 140 pc; Elias 1978). The primary and secondary stars have 
spectral types of M0 and M3.5, respectively, with a projected separation, a sep , of 0.242" (34 
AU; Hartigan & Kenyon 2003). The dynamical mass of the binary is calculated to be 0.78 ±0.25 
M (Tamazian et al. 2002). Based on the shape of the spectral energy distribution (SED), 
the binary is classified between Class I and II (Andrews & Williams 2005). 

Gledhill & Scarrott (1989) suggested a circumbinary disk with a diameter of ~ 5" (700 
AU). The disk shows a linear polarization of ~ 10 % in the i?-band, which is higher than that 
of its associated molecular cloud (~ 2 %; Vrba et al. 1985). However, Krist et al. (1998) could 
not confirm a disk beyond 60 AU from the binary using HST / Wide Field Planetary Camera 2 
(WFPC2) observations at the V-, R-, and J-band wavelengths, although they did not subtract 
the point-spread functions (PSFs) of the binary. The sub-millimeter and millimeter continuum 
fluxes of the binary are extremely low in comparison with those of other T Tauri stars (Andrews 
& Williams 2005; Beckwith et al. 1990; Dutrey et al. 1996). The gas and dust mass around 
the binary is calculated to be 0.002 M Q . 

FS Tau is actually a hierarchical triple system. The third component, called as FS 
Tau B or Haro 6-5 B, is located about 20" (2800 AU) west of the FS Tau binary. However, 
high spatial resolution imaging observations revealed no directly visible star in the Haro 6-5 
B system (Krist et al. 1998). Haro 6-5 B itself is obscured by its circumstellar disk, which is 
detected in millimeter wavelengths (Dutrey et al. 1996). The companionship of Haro 6-5 B to 
FS Tau A is not proven by common proper motion, spectroscopy, nor photometry. Instead, the 
companionship is indicated by polarimetry. The FS Tau A binary is surrounded by a reflection 
nebula up to 1400 AU in radius, which shows a centro-symmetric polarization pattern centered 
on the binary, particularly on the south side (Gledhill & Scarrott 1989). Krist et al. (1998) 
suggested that the north side is illuminated not only by the FS Tau binary but also by Haro 6-5 
B. H a imaging observations indicated a cavity wall in the north side, excited by the blue-shifted 
outflow from Haro 6-5 B (Eisloffel & Mundt 1998). The reflection nebula and the cavity wall 
appear to be truncated by the straight edge of a dark cloud ~ 700 AU north of the binary (Fig. 
2; see also Fig. 3 of Krist et al. 1998). 

The FS Tau binary has not been known to drive outflows. Optical emission lines, such 
as [S II] and H Q , showed no evidence of high- velocity jets from the binary (Woitas et al. 2002). 
Although 13 CO emission was observed along the southeast direction from the binary, it is 
unclear whether this emission arises from a diffuse wind (Dutrey et al. 1996). 

In this paper, we present a subarcsecond NIR image of the FS Tau binary taken with the 
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Coronagraphic Imager with Adaptive Optics (CIAO) on the Subaru Telescope. Combined with 
the HST optical imaging and polarimetric methods, we discovered complex features associated 
with the binary. Our goal was to dissect the circumbinary disk, whose southeast side is bright 
in NIR wavelengths, in contrast with the northwest side, which is bright in the optical. We 
interpret the difference in brightness as being caused by the forward scattering of dust on the 
southeast side and a contribution to the radiation from Haro 6-5 B on the northwest side. 
Observations and data reduction are described in §2. Our results and some discussion are 
presented in §3 and §4, respectively. Finally, our conclusions are presented in §5. 

2. Observations and Data Reduction 

2.1. CIAO 

Near-infrared (i/-band; 1.6 fim ) coronagraphic observations of the FS Tau binary were 
carried out on 2007 December 17 with the Subaru/CIAO (Tamura et al. 2000) system. A 0"8 
mask blocked the light from both the primary and secondary stars (a S cp ~ 0('24). This mask 
had ~ 2% transmission (Murakawa et al. 2003). A pupil Lyot stop reduced the diffracted 
light. The field of view was 21. "83 x 21."83, with a pixel scale of 0."02132±0."00003 pixel" 1 and 
a position angle (P. A.) on the detector of 1?58±0?07. The pixel scale and orientation were 
measured with observations of the Trapezium cluster (Simon et al. 1999). One may doubt 
whether the astrometrical calibration is still valid, since the data of Simon et al. (1999) were 
taken in 1996. The pixel scale and the P.A. derived here are consistent with those derived by 
Itoh et al. (2005), whose data were taken in 2003. Thus we believe that these parameters 
are correctly derived. During the observations, the natural seeing size varied between 0'.'9 and 
l'.'A in the if-band. The average full width at half-maximum (FWHM) of the PSFs of our 
observations was 0."4, using AO. We obtained 36 object frames in the if-band. The integration 
time was 10 s x 3 co-adds for each object frame. After obtaining the first 24 frames, we 
performed dithering to remove hot and bad pixels. We did not obtain science - sky - science 
jittering as it is the standard for infrared observations of extended sources. This is due to very 
low observing efficiency. For the coronagraphic observations with the AO system, the setup 
procedures are very complicated, as about 10 minutes was required for every pointing. After 
the FS Tau observations, we also obtained 25 frames of a PSF reference star, SAO 76648, in the 
H-band using the 0"8 mask. It is located about 140' east to FS Tau. The //-band magnitude 
of the PSF reference star was only 0.5 mag brighter than that of FS Tau. On the other hand, 
the i?-band magnitude for the wavefront sensing is 1.5 mag brighter than that of FS Tau. We 
adjusted it to match FS Tau using neutral density (ND) filters in the AO. This was necessary 
to make the PSF of the reference star similar to that of FS Tau. The integration time was 10 
s x 3 co-adds for each reference star frame. FS 123 was observed as a photometric calibrator 
(Hawarden et al. 2001). Twilight-flats and dark frames were taken at the end of the night. 
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The object frames were calibrated using the Image Reduction and Analysis Facility 
(IRAF). First, a dark frame was subtracted from the object frames. Then the frames were 
divided by the normalized twilight-fiat to remove pixel-to-pixel variation in sensitivity. Hot 
and bad pixels were removed from the divided frames using the COSMICRAYS task. Then, an 
averaged count in the sky region was subtracted. We judged that neither halo of the PSF nor 
circumstellar structures extended at 6" away from the central star. The count was measured 
along a concentric circle (radius of ~ 6" with a width of ~ 1" ) on the FS Tau binary, avoiding 
the diffraction artifacts of the spider. The same reduction procedures were applied to the 
reference star frames. 

We subtracted the reference star frames from the object frames to detect faint structures 
buried in the halo of the central binary. The peak positions of the FS Tau binary and reference 
star were measured with the IMEXAMINE task on all frames. For FS Tau, the peak positions 
were able to be measured for 21 frames. For the rest frames, we did not identify the core 
of the PSF. The secondary star (FS Tau Ab) was not resolved in the object frames because 
the averaged FWHM was about 0''4. We assumed that the measured positions of the PSF 
peak corresponded to the positions of the primary. The object frames were shifted so that the 
positions of the primary were centered on the image. The reference star frames were duplicated 
after dark-subtraction, flat-fielding, hot and bad pixel rejection, and sky-subtraction. Then 
each reference star frame was shifted to adjust the peak positions of the PSFs between the 
reference star and each component of the binary. We calculated the positions of the secondary 
from the separation and P. A. in the H ST /ACS image (see §3.2). As the binary has only a 
separation of 0."23 (32 AU) and its total mass is 0.78 M this yields an orbital period of about 
210 years. Hence, the position angle of the binary could have changed up to 5° within the epoch 
difference of 3 years between the ACS and CIAO observations. It produces 1 pixel uncertainty 
(0."23 x tan 5° ) in the position of the secondary star. The intensities of the reference star frame 
were scaled to adjust the flux ratio of the binary. Because the if -band flux ratio of the binary 
is unknown, we presumed that it was the same as the fT-band flux ratio (7.0 ±0.1; White & 
Ghez 2001). One may consider that the if -band flux ratio can be calculated from the spectral 
types of Aa and Ab and the fT-band flux ratio. However, since FS Tau has near-infrared excess 
(J = 10.705 ±0.027 mag, H = 9.244 ±0.026 mag, K = 8. 178 ±0.017 mag from 2MASS) and we 
cannot identify whether the excess arises from FS Tau Aa and/or Ab, we did not determine 
the flux ratio of Aa and Ab in the if -band. We are convinced from the final PSF-subtracted 
image that the assumed flux ratio is appropriate. Then the 12 reference star frames, except 
the frames with the asymmetric PSFs, were combined. The combined reference star frames 
were renormalized to match the intensities between the reference star frames and the object 
frames in the 0'.'5x0''5 region immediately outside the mask southwest of the binary, where 
the effect of the spider was negligible. The reference star frames were subtracted from each 
object frame. Finally, the 21 PSF-subtracted object frames, except the frames. We checked 
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the photometric accuracy under the unstable seeing condition. Any object except for FS Tau 
itself did not appear in the object frames. FS Tau in the object frames is not appropriate for 
evaluating the photometric accuracy, because the core of the PSF was under the mask. Instead, 
we evaluated it with the i?-band flux of FS Tau, which was monitored during the integrations 
by the wavefront sensor in the AO system. The measured fluxes were stable with its standard 
deviation of 0.04 mag. We consider that our observations have photometric accuracy of 0.04 
mag. 

2.2. HST Polarimetric Data 

Polarimetric observations were carried out on 2004 August 21 with the HST Advanced 
Camera for Surveys (ACS). These observations were proposed by D. Hines (HST program 
GO 10178). The central wavelength and the width of the F606W filter (broad V^-band) are 
5907 A and 2342 A , respectively. The field of view was 202" x202" , with a pixel scale of 
0"050±0"005 pixel -1 (The ACS manual). The position angles of the images were recorded in 
the FITS keyword of the image, ORIENTAT. The uncertainty in the position angle was 0?01 
- 0.°03 in the worst cases (Koekemoer et al. 2007). The exposure time of each polarized image 
(0° , 60° , and 120° ) was 536 s. The averaged FWHM of the observations was 0"09. 

We downloaded the polarimetric data with the basic calibration from the HST MAST 
web site. Each frame was shifted so that the positions of the unsaturated FS Tau primary were 
centered on the image. The Stokes parameters, /, Q, and U were computed using 

J = ^(F + F 6 o + F 120 ) (1) 

Q = ^(2F -F 60 -F 120 ) (2) 

U = ^=(F m -F 120 ), (3) 

where F , F 60 , and F U o are the counts of each polarized image (Biretta et al. 2004). The degree 
of linear polarization P and the polarization angle 9 were computed from 

P = ^+E (4) 

9 = -arctan(^). (5) 
2 Q 

The uncertainty in the degree of polarization was ~ 0.3%, and the uncertainty in polarization 
angle was < 2.°0. 

We subtracted the PSFs of the primary and secondary stars in the Stokes / image from 
those created by the Tiny Tim program (model PSFs; Krist 2004) to detect circumbinary 
structures using the same data reduction as CIAO. The intensities of the model PSFs were 
scaled to adjust the ^/-band flux ratio of the binary (20.8; Krist et al. 1998). We combined 
the model PSF frames, and then the combined frame was renormalized to match the intensities 
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between the model PSF frames and the object frames of the 0"5xO"5 southwest region of the 
binary, immediately outside the "artifact region" (see §3.2). 

3. Results 

3.1. H-band Coronagraphic Image 

The //-band coronagraphic image is shown in Fig. 1. Despite bad seeing conditions 
an image with 0''4 FWHM could be obtained. The FS Tau binary exhibited four remarkable 
structures: 

1. A bright portion extending 630 AU southeast from the binary ("1" in Fig. 1). 

2. Two arm-like structures on the western side ("2" in Fig. 1). 

3. Two cavities at the northeast and southwest at P. A. ~20° and ~250° , respectively ("3" 
in Fig. 1). 

4. A faint (>17.3 mag arcsec -2 ) arc-like structure ~900 AU west of the binary ("4" in Fig. 
1)- 

We measured the surface brightness of the extended features with 0'/2 x 0'/2 square 
apertures. The southeast bright structure has 15.2±0.3 mag arcsec -2 at 220 AU from the 
central star and 16.8±0.2 mag arcsec -2 at 530 AU. If the PSF subtraction of the reference star 
were too strong or too weak, the residual would be seen as a circle or streak. We confirmed 
that the detected structures are real by making a pseudo-binary frame from the first half of the 
reference star frames. The projected separation and flux ratio of the pseudobinary were the 
same as those of the FS Tau binary (§2.1). We subtracted the PSFs using the same procedure 
as for FS Tau. In the processed image, we identified subtraction residuals within l."3 — 1."7 from 
the central star (Fig. 1). We called this region the "artifact region." Beyond this region, the 
diffraction pattern of the spider was the only artifact. We determined the limiting magnitude 
of the observations by measuring standard deviation of the sky regions which was set far from 
the central star and the nebulosity. We found that the limiting magnitude with S/N=3 was 
17.3 mag arcsec -2 . 

3.2. F606W-band Stokes I Image 

Figure 2 represents the F606W-band Stokes / images of the FS Tau/Haro 6-5 B field 
and the closeup FS Tau binary. The primary and secondary stars of the FS Tau binary were 
resolved in these images. The projected separation and P.A. were 0."23±0."01 and 104.°69±6.°0, 
respectively. The F606W-band image reveals the western arms ("a" in Fig. 2), the northeast 
and southwest cavities ("b" in Fig. 2), and the arc- like structure ("c" in Fig. 2) in common 
with the if-band coronagraphic image. Also, we discovered two eastern arms ("d" in Fig. 2) 
in the F606W-band image. One is sharp, but the other is not. These western and eastern arms 
cover the southwest and northeast cavities with opening angles of ~45° and ~40° , respectively. 
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A local peak was identified at l'.'O (140 AU) from the binary in the F606W-band images 
before and after the PSF-subtraction ("e" in Fig. 2). The FWHM was ~0"8 (110 AU). We 
identified the local peak as an extended dusty structure, not a stellar object (§4.2.2). Its 
F606W-band magnitude was measured to be 16.0±0.5 mag arcsec -2 with a 0."4 radius aperture; 
however, it is located in the artifact region and on the diffraction pattern of the spider. 

To evaluate the PSF-subtraction in the HST image, a pseudo-PSF binary was created 
using a dwarf, 2MASS J07040548-0350493 (HST program GO 9588 proposed by H. Bond). 
The flux ratio of the pseudo-binary was set to that of the FS Tau binary (20.8; Krist et al. 
1998). The projected separation of the pseudo-binary had the same value as the FS Tau binary 
in the F606W-band image. We subtracted the PSFs created by the Tiny Tim from the pseudo- 
binary. The PSF-subtraction produced artifacts within about 1% from the central star (Fig. 
2). The limiting magnitude of the observations was calculated from standard deviation of the 
sky regions, which was set far from the central star and the nebulosity. We derived that the 
limiting magnitude with S/N=3 was 22.1 mag arcsec -2 . 

3.3. F606W-band Polarimetric Image 

Figure 3 represents the F606W-band polarized intensity image and if-band corona- 
graphic image, together with polarization vectors. To clarify the sources of the polarizations 
around the FS Tau binary, the polarization angles were investigated. We defined the devia- 
tion from the centro-symmetric polarization angle centered on the binary and Haro 6-5 B as 
#centroA = \Q — P-A.a±|, and C cntroB = \6 — P.A.bjJ. These angles represent the observed polar- 
ization angles; P.A.aj. and P.A. B j_ are perpendicular to the P.A. with respect to the FS Tau 
primary and Haro 6-5 B, respectively, at each position. The red-painted polarization vectors 
in Figure 3 represent ^centroA < 30° and ^centroB > 30° ; the blue-painted polarization vectors 
#centroA < 30° and 6> ccntroB < 30° ; the green-painted polarization vectors # C entroA > 30° and 
#centroB < 30° ; the yellow-painted polarization vectors # C entroA > 30° and C entroB > 30° . If the 
FS Tau binary were the only illuminating source in the FS Tau/Haro 6-5 B field, the angles 
#centroA would fall below 30° at all points. This means that the blue-painted polarization vec- 
tors would be present between the P.A. = 57° to 135° and 250° to 303° if the binary had a 
geometrically flat circumbinary disk with a major axis along a P.A. = 30° and an inclination 
by 35° (§4.1.1). The red-painted polarization vectors would be present at the other P.A.s. 
Their degrees of polarization depend on the scattering angles. 

The FS Tau binary system has two polarimetric characteristics: small # C entroA features 
at the local peak (# ce ntroA = 5° ±4°), the southwest arm (7° ±5°), and the western cavity 
(13° ± 11°); and large # ce ntroA features at the portion of the southeast that is bright in the 
if -band (30° ±22° ) and at the north side at a distance from the binary, r, of r > 420 AU 
(53°±17° ). Causes of these large deviations in the polarization angles will be discussed in the 
next section. The polarized image does not show parallel polarization vectors or null points near 
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the central binary both of which are indicative of a polarization disk. The averaged F606W- 
band polarization was P = 11 ±5% in r < 560 AU, which is consistent with the result of Gledhill 
& Scarrott (10% at the i?-band). 

4. Discussion 

4-1. Circumbinary Disk 

We consider that the emission nebulae are mostly attributed to the circumbinary disk. 
4-1.1. Optical and NIR Imaging /Polarimetric Characteristics of the Disk 

Figure 4 shows the if -band radial profile along the bright southeast region. The if -band 
brightness was 15.2±0.3 mag arcsec -2 at r = 220 AU and 16.8±0.2 mag arcsec -2 at r = 530 
AU. We fit a single power-law of p -1-9 1 to the southeast region. Its slope was shallower 
than that of circumstellar disks around other classical T Tauri stars, such as FN Tau (-2.5; 
Kudo et al. 2008), TW Hya (-2.6; Weinberger et al. 2002), and GM Aur (-3.5; Schneider et 
al. 2003), in the ii-band. The surface brightness of a geometrically flat and optically thick 
structure fits a single power-law of r -3 , whereas the surface brightness of a flared structure with 
well-mixed gas and dust fits a power law of r -2 (Whitney & Hartmann 1992). We argue that 
the southeast region corresponds to a part of the circumbinary disk and is more flared than the 
other T Tauri disks. One may consider that the structure is a flattened envelope rather than 
a disk. But the central binary is not heavily reddened. The combined near-infrared colors are 
J — H = 1.46 ± 0.04 mag and H — K = 1.07 ± 0.03 mag, corresponding to the intrinsic colors of 
classical T Tauri stars (Meyer et al. 1997) with a moderate extinction (Ay ~ 6 mag). Based on 
these colors, we consider that the structure is a circumbinary disk. The ii-band radial profile 
of the northwest region of the binary are also shown in Figure 4. It is fainter than the southeast 
region at all distances. We fit a single power-law of r -L5±0 - 3 to the northwest region. Its slope 
seems shallower than that of the southeast region, but its faintness prevents us from further 
discussion. 

Figure 5 represents the F606W-band azimuthal profiles of the degrees of polarization 
around the FS Tau binary. The degrees of polarization were derived from a 0"25x0"25 aper- 
ture. We determined that the outer radius of the disk is consistent with that of the south- 
east region (4'.'5 = 630 AU), which is bright in the ii-band. The observed polarizations of 
#centroA > 30° (i.e., the green- and yellow-painted vectors in Fig. 3) were not applied in or- 
der to exclude the polarizations caused by multiple scattering events from dust grains in the 
southeast region (§4.1.4). In addition, the polarizations in the r > 3" north (P.A. = 280° to 
55° ) were masked due to the contamination by light from the cavity wall associated with Haro 
6-5 B (§4.1.2). The maximum degrees of polarization, P max , were observed in the northeast 
and southwest sides of the disk, and the minima, P m i n , were observed in the southeast and 
northwest sides. 
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We constructed a simple circumbinary disk model to investigate the polarization by 
dust grains in the FS Tau binary system. The FS Tau binary has a prominent smooth emission 
feature at 10 /im (Furlan et al. 2006), which indicates sub-micron-sized amorphous silicate 
dusts in the disk surface layer. Such dust grains show Rayleigh-like scattering. The GG Tau 
binary also shows a strong 10 /im silicate emission feature similar to that of the FS Tau binary 
(Furlan et al. 2006). The circumbinary disk around GG Tau displays a high polarization 
amplitude (~50%), with a centro-symmetric pattern in the NIR wavelengths (Silber et al. 
2000), which is well reproduced by Rayleigh-like scattering (a <C A) from sub-micron dust grains 
in the disk. Therefore, we assume that "single" Rayleigh-like scattering occurs in the FS Tau 
circumbinary disk. Using small dust grains with a power-law size distribution of n(a) oc a -3 5 
between 0.005 — 0.25/im (hereafter MRN dust; Mathis et al. 1977), the degree of polarization, 
P, depends on the scattering angle, 6* sca t, as below: 



where I\\ and I± are the intensities of scattered light in polarization modes parallel and per- 
pendicular to the scattering plane, respectively (van de Hulst 1957). The symbol i is the disk 
inclination, cf> is the flared angle of the disk, and x, y are the coordinates deprojected along the 
major and minor disk axes, respectively, with their origin at the center of the primary star (Fig. 
5), i.e., +0 for the near side of the disk along our line of sight and — for the far side. Three 
parameters, i, <p, and the P.A. of the disk major axis, £, are free parameters. These parameters 
were determined from a x 2 minimization method at r = 1" intervals. 

The disk model shows a centro-symmetric polarization pattern around the binary. The 
-Pmax occurs at two points, each at a 90° scattering angle (i.e., disk major axis). The disk 
minor axis corresponds to the positions of the P m i n , where the scattering angles are far from 
90° . If the disk is flared, P m i Q occurs at the near side of the disk. Notice that a spherical dust 
distribution, such an envelope, does not produce a sinusoidal pattern; in that case, the degree 
of polarization is constant as a function of the scattering angle (or P.A.). Because the observed 
P max occurs at P.A. = 30° and 210° at r = 2 - 3" (280 - 420 AU; Fig. 5 and Table 1), we 
determined the disk major axis at P.A. = 30° and a disk inclination of 35° . The flat disk 
model (<f) = 0°) agree well with the observed polarization at r = 2 — 3" . The polarization is at a 
minimum at P.A. = 120° and 300° . Assuming the southeast side is brighter than the northwest 
side in the if -band due to the forward scattering of dust, the southeast side is located closer 
to us (6> scat ~ 55°), and the northwest side is oriented away from us (9 scaut ~ 125°). However, 
the radial profile of the surface brightness of the disk is not consistent with the flat-disk model. 
The x 2 minimization method shows that the southeast side of the disk has a lower degree of 
polarization than the northwest side at r = 4 — 4.5" (560 — 630 AU). The flared disk model 




I± ~ I\\ _ 1 - cos 2 6> scat _ 
I± + I\\ 1 + cos 2 6» scat ' 
ytan(z±0) \ 



(6) 



(7) 
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(<f) = 5°) possibly explains the observed polarization at r = 4 — 4.5" , where the disk is inclined 
by 30° , with its southeast side toward us. 

The P.A. of the disk major axis and the inclination of the disk are consistent with 
those of the binary orbit, whose major axis is along the P.A. = 18°±7° and an inclination of 
20°±5° (Tamazian et al. 2002). The y 2 minimization method also shows (£, i, <ft) = (15° , 
40° , 0° ) at 1 - 2" and (35° to 40° , 30° to 35° ,0° to 5° ) at 3 - 4" . 

In conclusion, the circumbinary disk is characterized by the following geometry: the 
radius of the disk is 630 AU (4."5), the disk major axis is along the P.A. = 15° to 40° , and the 
disk is inclined by 30° to 40° , where the southeast side (P.A. = 105° to 130° ) corresponds to 
the side nearest to us. The outer (>560 AU) portion of the disk is possibly more flared than 
its inner portion. 

4-1-2. Southeast to Northwest Brightness Asymmetry 

The circumbinary disk shows an unambiguous southeast-northwest asymmetry. The 
brightness ratio of the southeast and northwest sides at r = 280 — 630 AU is 1.6 ± 0.4 in the 
if -band. We suggest that anisotropic scattering of the dust causes the asymmetric brightness 
of the disk and that the southeast side is brighter due to forward scattering of the dust. This 
indicates that the southeast side corresponds to the near side. Based on this hypothesis, the 
brightness ratio (B nc3iT / 5f ar ) of the disk is described as: 



where g is the asymmetry parameter of the Henyey-Greenstein scattering phase function 
(Henyey & Greenstein 1941), #Nscat is the averaged scattering angle of the near side, and #Fscat 
is that of the far side (Eq. 7). For the disk parameters in §4.1.1 (r = 630 AU, £ = 15° to 
40° , i = 30° to 40° , and 0=0° to 5° ), the ii-band brightness ratio was fit with g ~ 0.2. 
This value is consistent with that of the MRN dust (g ~ 0.1) shown by Kim et al. (1994), who 
assumed single scattering by the dust in a diffuse interstellar cloud (i?v = 3.1). For the FS Tau 
circumbinary disk, such small dust grains probably produce single-scattered light, as the disk 
is optically thin in the if -band (§4.1.4). 
4-1-3. Northwest Side of the Binary 

The northwest side of the circumbinary disk are bright in the optical wavelengths but 
faint in the near-infrared wavelengths (Fig. 6). We measured its color. Within l'/5 from 
the central star, the PSF subtractions are not perfect either in the optical image nor in the 
near-infrared image. Beyond 4" from the central star the disk is very faint in the near- 
infrared wavelengths. We measured the surface brightness of the disk in the region between 
l'/5 and 4'/0 from the central star and the P.A. between 280° and 355° . The F606W-if 
color in this region was derived to be between 1.7 mag and 3.0 mag. We also derived the 
intrinsic colors, (F606W— H) int , of the photosphere of the binary and Haro 6-5 B as 2.44 ±0.02 




3/2 



(8) 
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mag and 2.00 ± 0.02 mag, respectively, using the relationship (F606W— H)i at — (V — H) int — 
(0.477 ± 0.013) {V-R) int + (0.001 ± 0.005) (V-R)f nt (Sirianni et al. 2005). The intrinsic colors, 
(V — -R)i n t and (V — H) int , were derived from those of main sequence stars (Ducati et al. 2001) 
whose spectral types were the same as the FS Tau primary (M0; Hartigan & Kenyon 2003) 
and Haro 6-5 B (K5; White & Hillenbrand 2004). Comparing the disk color with the colors of 
the point sources, we cannot identify whether the disk is illuminated by the FS Tau binary or 
Haro 6-5 B. 

On the other hand, the polarization map clearly indicates the contribution of Haro 6-5 
B. We claim that the scattered light and optical emission lines from Haro 6-5 B change the 
degree of polarization and polarization angle on the north side at r > 3" . Several observations 
identified a cavity wall in the north side, which is created by the blue-shifted outflow from Haro 
6-5 B (Gledhill & Scarrott 1989; Eisloffel & Mundt 1998). Gledhill & Scarrott (1989) reported 
that the cavity wall has non-centro-symmetric polarization patterns with respect to the FS Tau 
binary, bluer colors than those of the south side of the binary, and an H a emission line. We 
confirm the radiation from Haro 6-5 B at r > 3" north (P. A. = 280° to 55° ) in the F606W-band 
polarimetric image. The green-painted polarization vectors are predominant (56% of the total) 
in the r > 3" north side (Fig. 3). The averaged # C entroA is 53° ±17°, whereas the averaged 
#centroB is 11° ±8°. The averaged P is 16 ±7%. These ^centroB and P are comparable to those 
of the cone-shaped reflection nebula (hereafter Rl nebula) that extends 8" northeast from its 
illuminating source Haro 6-5 B. In addition, the radiation from Haro 6-5 B is a possible cause 
of the blue-painted polarization vectors at r > 3" north. The averaged ^centroA is slightly larger 
than the averaged # cent roB (20° ± 7° vs. 14° ±9°). We argue that the cavity wall blocks the 
north side of the circumbinary disk at r > 3" . On the other hand, for the north side at r < 3" , 
the averaged ccntTOJ ± is 13° ± 13° and the averaged ccntroB is 44° ±24°, indicating the scattered 
light from the FS Tau binary. 
4- 1-4- Southeast Side of the Disk 

The southeast side of the circumbinary disk is bright in the near-infrared wavelengths 
but faint in the optical wavelengths (Fig. 6). We measured its color in the region between l."5 
and 4."5 from the central star and the P.A. between 93° and 175° . The F606W— H color of this 
region of the disk is derived to be 4.2±0.2 mag, which is redder than the intrinsic color of the 
FS Tau binary (Fig. 6). Figure 3 displays a weak centro-symmetric (# C entro ~ 30°) polarization 
pattern around the southeast portion, despite the high signal to noise ratio (> 6). We consider 
that the southeast portion corresponds to a part of the circumbinary disk. 

Given the (F606W-#) int of the binary (2.44 ±0.02 mag) and the F606W-# color of the 
southeast portion (4.2±0.2 mag), the color excess of the southeast portion is -E(F606W— H) = 
1.8 ±0.2 mag. We assume that extinction by dust grains around the binary reddens the south- 
east portion. This extinction, A\, was derived as 2.3±0.3 mag in the F606W-band and 0.5±0.1 
mag at the if-band. We estimated an optical depth, t\, of 2.1 ±0.3 in the F606W-band and 
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0.5 ±0.1 at the if-band, using the relationship A\ = 1.086rx (e.g., Stark et al. 2006). Thus, the 
southeast portion is optically thick in the F606W-band, and some of the scattered light comes 
from the southeast portion, where it expected multiple scattering events from dust grains. 
Further circular polarization observations are needed to perceive the multiple scattered light 
events. 

4.2. Outflow 

No optical jet was detected within a few arcseconds of the FS Tau binary (Woitas et al. 
2002). Although 13 CO emission was observed along the southeast side, it is unknown whether 
this emission arises from the diffuse wind from the binary (Dutrey et al. 1996). The circumbi- 
nary disk model (§4.1.1) predicts a red-shifted outflow driven by the primary or secondary stars 
along the southeast direction (P. A. = 105° to 135° ), given that the circumstellar disks are 
coplanar with the circumbinary disk, and the outflow axis is perpendicular to the circumstellar 
disk plane. 

In contrast, both the if -band and F606W-band images show evidence of other outflows: 
the western and eastern arms + cavity systems, and the arc-like structure (Figures 1 and 2). 
The arc-like structure located at r ~ 900 AU west is seen faintly in the if -band image, but is 
clearly detected at above 10<r in the F606W-band image. This structure appears to connect 
the western arms or to have an appearance similar to bow shocks excited by Herbig-Haro flows. 
The two western arms encompass the southwest cavity with an opening angle of ~45° . The 
eastern arms + cavity system, only detected in the F606W-band image, plausibly mixes with 
the cavity wall associated with Haro 6-5 B (§4.1.3). We argue that these arms + cavity systems 
are similar to the features created by bipolar outflows in several T Tauri systems (e.g., Haro 
6-5 B: Eisloffel & Mundt 1998, HL Tau; Lucas et al. 2004). We measured surface brightness of 
the cavities. The regions of the cavities were defined as circular areas with 3. "5 radius centered 
at 12"4 with the P.A. of 43?5 from the central star and centered at 9"9 with P.A=245?7 for 
the northeast cavity and the southwest cavity, respectively. The (b) marks of the right figure of 
Fig. 2 correspond to the central positions of the cavity regions. The southwest and northeast 
cavities have comparable F606W-band brightnesses of 21.2 ±0.2 mag arcsec -2 , whereas their H- 
band brightness are lower than the limiting magnitude (17.3 mag arcsec -2 ). Their F606W— H 
colors are < 3.9 mag i.e., tf606W < 1-9 and th < 0.4 (see §4.1.4), assuming that the cavities are 
reddened by dust extinction. The outflow from the binary would sweep up materials in the 
cavities. The polarization of the northeast cavity is 18 ±5%, and that of the southwest cavity 
is 19 ±6%. These polarizations are relatively higher than the other structures located within 
560 AU of the binary (11%). The southwest cavity shows 6* centro A = 13° ± 11°, indicating the 
efficiently scattered light from the binary. The northeast cavity has ^centroA = 22° ± 15° and 
#centroB = 22° ± 10°. The northeast cavity is illuminated by mixed radiation from the binary 
and Haro 6-5 B, although the dereddened color of the mixed radiation is unclear. We claim 
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that the outflow emanates to the P.A. = 220° to 250° , which is equal to the axisymmetric 
orientation of the western arms. This direction is not consistent with that calculated by our 
disk model (P.A. = 105° to 135° or 285° to 315° ). These differences are possibly produced 
by the mechanisms discussed in the following sections. 
4-2.1. Disk Misalignment 

One possible mechanism is that the primary /secondary circumstellar disks are not copla- 
nar with the circumbinary disk. For T Tauri binaries separated by several AU (e.g., DQ Tau; 
Mathieu et al. 1997, UZ Tau E; Prato et al. 2002) and T Tauri binaries with moderate separa- 
tion (e.g., GG Tau; Tamazian et al. 2002), the circumbinary disks are coplanar with each binary 
orbital plane. In contrast, for widely separated T Tauri binary systems, i^-band polarimetric 
observations suggest that circumstellar disks are not perfectly coplanar but differ slightly from 
each other (~ 20°; Jensen et al. 2004). 

Bate et al. (2000) demonstrated that during about 2000 binary orbital periods circum- 
stellar disks are roughly aligned with a binary orbital plane by tidal torques. This timescale is 
estimated to be 5.5 x 10 5 yrs for the FS Tau binary, given its separation of 39 AU and its total 
binary mass of 0.78 M & (Tamazian et al. 2002). We suggest, based on a comparison of the 
alignment timescale with the stellar age based on the shape of the SED (10 5 — 10 6 yrs; Andrews 
& Williams 2005), that the circumbinary disk is still misaligned with the circumstellar disks. 
4-2.2. Gravitational Instability and Inhomogeneous Distributions of the Disk Materials 

Another possibility is that the western and eastern arms are formed by other mecha- 
nisms, such as gravitational disk instability or inhomogeneous distributions of the circumbinary 
materials. 

A circumstellar disk associated with AB Aur has spiral arms that are expected to be 
formed by the gravitational instability of the disk (Fukagawa et al. 2004). Nelson et al. (1998) 
revealed that such a spiral arm appears in a circumstellar disk if Q < 2.0, where Q = c s Q/ (ttGH) 
is Toomre's Q-value, c s is the sound speed, Q is the angular velocity, and X is the surface density 
of the disk. However, the FS Tau circumbinary disk is less massive (0.002 M & ; Andrews & 
Williams 2005). Its Toomre's Q- value is much larger than 2.0, given the disk radius of 630 AU, 
the surface density distribution of X oc r _a5 (Kitamura et al. 2002), and c s = 190 m s^ 1 (T = 10 
K, for typical T Tauri disks). Thus, we conclude that gravitational disk instability does not 
occur in the circumbinary disk. 

The second possible mechanism is that a companion shapes the spiral arm-like features 
(e.g., Bate et al. 2003). However, high spatial resolution observations have detected no other 
companion around the FS Tau binary (Chen et al. 1990; Simon et al. 1992; Krist et al. 1998; 
Hartigan & Kenyon 2003; Connelley et al. 2009). We identified a local peak at r ~ 140 AU 
north in the F606W-band polarimetric images ("e" in Fig. 2 and "A" in Fig 3). This structure 
corresponds to the "small arc" detected by Krist et al. (1998). This local peak is extended (its 
FWHM of 110 AU) and shows P = 15% and # ccntro A = 5°. These indicate that the local peak 
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is an extended dusty structure associated with the binary, not a stellar companion. The lower 
limit of the local peak mass is estimated to be ~ 0.3 M e from its F606W-band brightness (16.0 
mag arcsec" 2 ). Such a low mass structure cannot create the disk instability. 

Finally, several studies have demonstrated arm-like structures as resulting from inhomo- 
geneous distributions of materials in a circumbinary disk (Giinther & Kley 2002; Ochi et al. 
2005). The observed local peak could be such an inhomogeneity. However, no inhomogeneous 
gas and dust structures around the binary are detected in the millimeter wavelength obser- 
vations due to insufficient spatial resolution (~ 3" ) and flux sensitivity (Dutrey et al. 1996; 
Andrews & Williams 2005). The Atacama Large Millimeter/sub millimeter Array (ALMA) will 
clarify whether the arms + cavity systems correspond to outflowing materials or to a part of 
the Keplerian circumbinary disk. 

5. Conclusions 

Using the Subaru/CIAO, we obtained an iJ-band coronagraphic image of the T Tauri 
binary system FS Tau. Combined with the H ST / ACS (F606W) polarimetric images, several 
features around the binary were identified. 

The circumbinary disk extends to 630 AU in radius. The disk is inclined by 30° to 
40° , such that the southeast side is nearest. The outer portion of the disk may be more flared 
than the inner portion. The disk also displays an unambiguous southeast-northwest brightness 
asymmetry. The southeast side of the disk is significantly red (F606W— H = 4.2 ±0.2 mag). 
It shows a weak centro-symmetric (# ccn troA ~ 30°) polarization pattern in the F606W-band. 
These indicate multiple scattering events in the optically thick southeast side at F606W-band 
wavelengths. The northeast side of the binary is about neutral in color (F606W— H = 1.7 mag or 
3.0 mag) compared with the FS Tau binary and Haro 6-5 B. We argue that the northwest side of 
the binary at r > 420 AU (3" ) is blocked by the cavity wall associated with Haro 6-5 B. This idea 
is supported by the observed polarization angles (6> centro A = 53° ±17° and ^ntr-oB = 11° ± 8°). 
The r < 420 AU northwest side has 6> ccntroA = 13° ± 13° and # ccntroB = 44° ± 24°, indicating 
scattered light from the FS Tau binary. 

The western cavity encompassed by the two arm structures is detected in both the 
F606W-band and if-band images. The eastern cavity, with additional arms, is seen in the 
F606W-band image, but only the cavity appears in the if-band. We argue that these arms + 
cavity systems are created by the bipolar outflow from the binary along the P.A. = 220° to 
250° . Another interpretation is that these structures correspond to inhomogeneous distri- 
butions of materials in the disk. Further high spatial and velocity resolution observations at 
sub-millimeter and millimeter wavelengths (e.g., ALMA) are required to clarify the mechanisms 
of these structures. 
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Fig. 1. The -ff-band coronagraphic image of the FS Tau binary system (left). The PSFs of the central 
binary were subtracted. "1" indicates the //-band bright portion of the circumbinary disk; "2" indi- 
cates the arm- like structures; "3" indicates the cavities; and "4" indicates the arc- like structure. The 
position of the FS Tau primary is indicated by the cross symbol in the occulting mask with a diame- 
ter of 0"8 (112 AU; a circular symbol). The binary is not resolved due to bad seeing conditions. The 
pseudo-binary was created using the PSF reference star SAO 76648 (right). The PSF-subtracted image 
shows some artifacts within 1"3 - 1"7 from the central star, represented as hatched ellipses in both im- 
ages. Beyond this region, the diffraction pattern of the spider is the only artifact, indicated by dotted 
lines in the both images. The field of view (FOV) is 17"8xl7"8. North is up and east is to the left. 
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Fig. 2. The F606W-band Stokes I images of the FS Tau/Haro 6-5 B field before subtracting the PSFs 
of the central binary {left) and a closeup of the binary after subtracting the PSFs (right), "a" rep- 
resents the western arms, "b" the cavities, "c" the arc-like structure, and "d" the eastern arms. The 
artifact region is shown by the hatched ellipse. The directions of the spider pattern are indicated by 
dotted lines, "e" represents the local peak in the logarithmic image of the FS Tau binary at lower 
left of the top panel. The positions of the FS Tau primary and secondary stars arc indicated by 
the big and small cross symbols, or "Aa" and "Ab," respectively. The cone-shaped reflection neb- 
ula Rl extends to 8" northeast from its illuminating source Haro 6-5 B. The dashed lines in the top 
panel indicate the straight edge of the dark i cloud suggested by Krist et al. (1998). The FOVs are 
36''3x49"4 in the top panel and 18''8x22"6 in the bottom panel. North is up and east is to the left. 
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Fig. 3. The F606W-band polarized intensity image {top) and il-band coronagraphic image (bottom; same 
as Fig. 1) with the F606W-band polarization vectors. The vectors are 0''25x0''25 binned polarization 
vectors and are delineated by intensities above 3cr in each polarized image. The vectors are also classified as 
follows: #controA < 30° and # CO ntroB > 30° (red), 8 ccntlo A < 30° and 6> CO ntroB < 30° (blue), C entroA > 30° and 
#controB < 30° (green), cen troA > 30° and ce ntroB > 30° (yellow). 9 ccntm A and 6> contro B are defined as the 
deviation from the centro-symmetric polarization patterns with respect to the FS Tau binary and Haro 6-5 
B, respectively (see text). The polarization patterns are broadly centro-symmetric around the circumbi- 
nary disk, resulting from illumination by the FS Tau binary. " A" represents the local peak in the disk. The 
polarization vector is non-centro-symmetric around the H-b&nd bright portion of the disk (" 1" ). The FOV 
is 19"0x30"6 in the top image and 13"8xl6"5 in the bottom image. North is up and east is to the left. 
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Fig. 4. The _ff-band surface brightness profile of the southeast (filled circle; P.A. = 93° to 175° ) 
and northwest (asterisk; P.A. = 280° to 355° ) sides of the binary. The brightnesses are aver- 
aged in a 0'.'2x0"2 aperture at r = 220 — 630 AU. The error bars represent standard deviations. The 
solid line represents a r^ 1 - 9 * - 1 power-law, and the dashed line represents a j- -1 ^ - 3 power-law. The 
limiting magnitude of the observations (17.3 mag arcsec -2 ) is shown by the horizontal dotted line. 
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r = 1-2 arcsec 



r = 2-3 arcsec 




Fig. 5. The azimuthal polarization profiles for various distances from the binary, starting from the top 
left (r = 1 — 2" ) and ending at the bottom right (r = 4 — 4.5" ). The observed polarizations with 
#centroA < 30° are applied to exclude the polarization due to multiple scattering events in the circumbinary 
disk. The observed polarizations in the r > 3" north (P.A. = 280° to 55° ) are not shown due to 
the contamination of the light from Haro 6-5 B. The solid curves represent the best match of the flat 
(<f) = 0°) circumbinary disk model. The dashed curves represent the best match of the flared (<fi = 5°) 
disk, where the southeast side corresponds to the near side. The dotted curves are the best match of the 
flared (4> = 5°) disks, where the northeast side is the near side. These x 2 are shown in Table 1. is 
the P.A. of the disk major axis, which is perpendicular to the plane of the bottom left panel, "i" is the 
disk inclination, and 4» is the flared angle. No data are shown at r < 1" due to the large uncertainty. 
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Fig. 6. The F606W— H color imago of the FS Tau binary. The hatched ellipse 

marks the artifact region caused by the PSF-subtraction in both the F606W-band 
{green) and ff-band images [red). The directions of the spider pattern arc indi- 

cated by dashed lines. The FOV is I6'.5xl7'.5. North is up and east is to the left. 
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Table 1. Best fit parameters for the circumbinary disk. §: the P. A. of the disk major axis; i: the disk inclination; <j>: the flared 
angle. 



Distance 




i 





x 2 


[arcsec] 


[deg] 


[deg] 


[deg] 




Flat Disk 


1 — z 


15 


40 


U 


i.i 


2-3 


30 


35 





2.1 


3-4 


40 


35 





1.7 


4-4.5 


35 


30 





0.9 




Flared Disk (SE) 






1-2 


15 


40 


5 


1.2 


2-3 


30 


35 


5 


2.3 


3-4 


35 


30 


5 


1.7 


4-4.5 


35 


30 


5 


0.8 


Flared Disk (NW)* 


1-2 


15 


40 


5 


1.5 


2-3 


35 


35 


5 


2.4 


3-4 


50 


35 


5 


2.0 


4-4.5 


45 


35 


5 


1.0 



* The near side of the flared disk corresponds to the southeast (SE) or northwest (NW) sides. 
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